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We discuss space-time chaos and scaling properties for classical non-Abelian gauge fields discretized on a spatial 
lattice. We emphasize that there is a "no go" for simulating the original continuum classical gauge fields over 
a long time span since there is a never ending dynamical cascading towards the ultraviolet. We note that the 
temporal chaotic properties of the original continuum gauge fields and the lattice gauge system have entirely 
different scaling properties thereby emphasizing that they are entirely different dynamical systems which have 
only very little in common. Considered as a statistical system in its own right the lattice gauge system in a 
situation where it has reached equilibrium comes closest to what could be termed a "continuum limit" in the 
limit of very small energies (weak non-linearities). We discuss the lattice system both in the limit for small 
energies and in the limit of high energies where we show that there is a saturation of the temporal chaos as a 
pure lattice artifact. Our discussion focuses not only on the temporal correlations but to a large extent also 
on the spatial correlations in the lattice system. We argue that various conclusions of physics have been based 
on monitoring the non-Abelian lattice system in regimes where the fields are correlated over few lattice units 
only. This is further evidenced by comparison with results for Abelian lattice gauge theory. How the real time 
simulations of the classical lattice gauge theory may reach contact with the real time evolution of (semi-classical 
aspects of) the quantum gauge theory (e.g. Q.C.D.) is left as an important question to be further examined. 



There are some indications - and it would be 
a beautiful principle if it was true - that we have 
the "approximative laws and regularities" which 
we know, because they are infrared stable against 
modifications of "short distance physics" in the ul- 
traviolet EL These "approximative laws" , which at 
present accessible scales to the best of our knowl- 
edge comprise a sector consisting of the Stan- 
dard Model of quantum Yang-Mills fields based 
on the group S{U 2 x U 3 ) ~ [7(1) <g> SU(2) <g> SU{3) 
and a sector consisting of the gravitational inter- 
actions, mav_ (thus) be robust and stable in the 
sense of Reffl, but they are full of unstable and 
chaotic solutions! Indeed, both the gravitational 
field (the Einstein equations) and the Yang-Mills 
fields exhibit dynamical chaostm for generic solu- 
tions (solutions without too much symmetry) in 
regions where (semi) classical treatments are justi- 
fied and non-linearities of the interactions are non- 
negligible. In the present contribution we shall 
report on some observations aau concerning the 



lattice implementation of space-time chaos of clas- 
sical Yang-Mills fields, but let us first note some 
few physical motivations for studying dynamical 
chaos in classical Yang-Mills fields: 

(1) Despite the word "chaos" at first suggests 
something structureless it is really a study of ways 
in which "structure" may be generated (from var- 
ious initial configurations) during the evolution of 
the governing equations of motion. It is of in- 
terest - but so far only little is known - to know 
how structures form and evolve (e.g. the formation 
and evolution of embedded topological structures) 
during the time evolution of the Yang-Mills equa- 
tions. In the semi-classical regimes of the Stan- 
dard Model, the quantum solutions build around 
the classical solutions (and it is not unlikely that 
quantum fluctuations on top of classical chaos only 
enhance chaos n). 



"The remark that the full quantum theory for a system 
with bounded configuration space (and a discrete spectrum 
of quantum states) has no chaos due to quasipcriodic evo- 
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(2) For example, the possibility of baryon 
non-conservation (via the famous ABJ-anomaly) 
within the electroweak theory has achieved quite 
some attention recently, and it is related to a de- 
tailed understanding of the dynamics of the elec- 
troweak fields (as the Universe cools down). It 
is quite natural to speculate about a relationship 
between dynamical chaos and an activity of for- 
mation and destruction of embedded topologically 
interesting field configurations (such a relationship 
is, for example, well known for the complex Lan- 
dau Ginzburg equation, see also e.g. discussion in 
Ref. H which contemplates the relevance of a con- 
cept of topological turbulence of the gauge fields). 
We note that field configurations for which the rate 
of production of baryon number 

B~ / d 3 x Tr(F F*) (1) 

vanishes will span a surface of co-dimension one. 
Thus, in fact, most field configurations (in the hot 
electroweak plasma) will contribute to the right 
hand side of equation (0). For a particular class 
of field configurations which contributes to the 
baryon non-conservation see also Rem. 

(3) Fast equilibration processes which take 
place in heavy-ion collisions!!!] are very likely con- 
nected to non-linear chaotic dynamics; an idea 
which in principle dates back e.g. to Fermi, Pasta 
and Ulam, Ref. 0. 

As a further motivation of the study of clas- 
sical Yang-Mills chaos, we should also note that 
not many non-pcrturbativc tools are available to 
study the time evolution of quantum Yang-Mills 
fields, so the study of classical Yang-Mills fields is 
a natural starting point for semi-classical under- 
standing of the dynamics. 

In order to facilitate a numerical study of 
spatio-temporal chaos (formation of space-time 
structure) in inhomogeneous Yang-Mills fields one 
has in some way to discretize the space-time con- 
tinuum on which the Yang-Mills fields are defined. 
This can be achievecLin a way which breaks gauge 
invariance (see e.g. E-j) or in a gaugeJjavariant way 
which is called lattice gauge theorylia. 

lution of the wave function is of little importance to us 
since the number of degrees of freedom involved here are so 
large that the time it will take to evolve quasi-periodically 
through the states will be substantially larger than the life- 
time of the Universe. 



We shall refer toem for more details in the dis- 
cussion which will follow. The study of chaos in 
toy-models where the gauge fields are spatially ho- 
mogeneous was initiated by Sergei Matinyan and 
George Savvidy 113. An important new qualita- 
tive dynamical feature comes into play, however, 
when one considers the spatially inhomogeneous 
classical Yang-Mills equations: There is a never 
ending cascading of the dynamical degrees of free- 
dom towards the ultraviolet, generated by the time 
evolution of the Yang-Mills equations! (In spite 
of this "ultraviolet catastrophe" , the solutions are 
well behaved in the sense that there are no "finite 
time blow up of singularities" ) . The non- linear 
self-coupling terms which open up the possibility 
for a chaotic behavior in the classical evolution 
lead in the non-homogeneous case to the infinite 
cascade of energy from the long wavelength modes 
towards the ultraviolet. (Note, there is a priori no 
concept of temperature and in the case of Abelian 
(electromagnetic) fields this cascading would not 
show up dynamically, unless one couples the fields 
to charged particles). This tendency of the mode 
frequencies cascading towards the ultraviolet will 
completely dominate the qualitative behavior of 
the classical Yang- Mills equations, and the "ultra- 
violet catastrophe" has for some time beem em- 
phasized by us (cf. e.g. discussion in Rem) as 
a major obstacle to simulate the classical contin- 
uum Yang-Mills fields in a numerical experiment 
over a long time span. |] There is no mechanism, 
within the classical equations, which prevents this 
never ending cascading of the modes towards the 
ultraviolet. Nature needs h, the Planck constant, 
as an ultraviolet regulator. Indeed, both Abelian 
and non- Abelian gauge fields are implemented as 
quantum theories in Nature. 

Here we shall discuss the possibility of using 
a lattice cutoff in a purely classical treatment to 
regularize the equations. There will still be a cas- 
cading of modes towards the ultraviolet, i.e. to- 
wards the lattice cut-off, and this ultraviolet cas- 
cade will still dominate the dynamical evolution 
of smooth initial field configurations. However, in 
a lattice formulation of the Yang-Mills fields on 
a large but finite lattice, the phase space is com- 
pact for any given energy and thus the system can 

b This is an obstacle which, in our opinion, has received 
insufficient attention in the various studies attempting at 
discussing and modeling chaotic properties of spatially in- 
homogeneous classical Yang-Mills fields. 
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reach an equilibrium state among the modes (a 
'thermodynamic equilibrium'). The lattice regu- 
larization of the theory opens up for the definition 
of dynamic and thermodynamic properties, which 
are not defined in the classical Yang-Mills field 
theory without regularization. It could e.g. be 
ergodic (modulo constraints) with respect to the 
Liouville measure, in which case it makes sense 
to talk about its micro-canonical distribution and 
approximating this by looking at 'typical' classi- 
cal trajectories. The fundamental assumption of 
thermodynamics asserts that on average the two 
approaches give the same result if we have a large 
system, and we may then naturally introduce cor- 
relation functions and possibly a correlation length 
£ (measured in lattice units) of the system. 

One hopes to define a continuum theory if, by 
judicious choice of the parameters in the system, 
one obtains a physical correlation length in the 
limit when the lattice constant goes to zero, i.e. 

£(a,E(a),...)xa I f , as a -> 0. (2) 

In equation (||) the correlation length £ in the lat- 
tice system is a function of lattice model parame- 
ters such as lattice spacing a, average energy den- 
sity E(a), etc. Condition (||) implies that the cor- 
relation length diverges when measured in lattice 
units and only if this is the case do we expect the 
lattice system to lose its memory of the underlying 
lattice structure. 

We shall restrict attention to the lattice gauge 
theory in 3+1 dimensions based on the gauge 
group SU(2). We consider a finite size 3 di- 
mensional hypercubic lattice having N 3 points 
where nearest neighbor points are separated by 
a distance a > 0. The phase space is a fibered 
space where the tangent manifold of the Lie-group 
SU(2) is assigned to each of the links, i 6 A, con- 
necting nearest neighbor lattice points: To each 
i € A is associated a link variable £7, € SU(2) as 
well as its canonical momentuir|] Pj e T Ui SU{2). 
A point in the entire phase space will be denoted 

X = {Ui,Pi} ieA e M = \\TSU(2) . (3) 

A 

A Hamiltonian is constructed so as to correspond 
to the continuum, classical Yang-Mills Hamilto- 
nian in the limit a — > 0. The construction of such 

c In principle Pi is a cotangent vector, but the Lie al- 
gebra inner product gives a natural identification of the 
cotangent space with the tangent space. 



a Hamiltonian is of course ambiguous in the sense 
that extra terms of order O(a) may be added to 
the Hamiltonian. Alluding to some sort of "univer- 
sality" (especially in the limit a — > 0), we expect 
that the precise choice of Hamiltonian is not so 
important in what follows. A Hamiltonian which 
is often employed to generate the time evolution of 
the orbit X(t) is the Kogut-Susskind Hamiltonian, 
which can be written in the following way |] : 

H(a,X^) = 

Here the last sum is over elementary plaquettes 
bounded by 4 links, and U a denotes the path- 
ordered product of the 4 gauge elements along the 
boundary of the plaquette □. The last term, the 
potential term, is automatically bounded and, for 
a given finite total energy, the same is the case for 
the first term, the kinetic term. Thus the phase 
space corresponding to a given energy-surface is 
compact. 

As is standard we shall discuss the temporal 
correlations (temporal chaos) of the lattice gauge 
system in terms of its spectrum of Lyapunov ex- 
ponents. The compactness of the phase space im- 
plies □ that the spectrum of Lyapunov exponents 
(which we overall will assume to be well defined 
quantities for the lattice system) is independent 
of the choice of norm on the_j3Dace of field con- 
figurations. In fact, it followsH'cl from the scaling 
properties of the equations of motion generated by 
the lattice Hamiltonian (||) , that in order to study 
the dependence of the maximal Lyapunov expo- 
nent with the energy density of the system, it is 
sufficient to consider the equations of motion for a 
fixed value of the lattice constant a, e.g. a = 1, as 
a function of energy density (cx energy/plaquette) 
and then rescale the results back afterwards. 

We have several different forms of lattice arti- 
facts in the lattice simulation of real-time dynami- 
cal behavior of the continuum classical Yang-Mills 
fields: 

(1) Lattice artifacts due to the compactness 
of the group. The magnetic term (the second 
term) in the Kogut-Susskind Hamiltonian (Q) is 

d For a derivation we refer to0 andFlfl from which we 
adapt our notation. For simplicity we omit the coupling 
constant factor 2/g 2 which anyway is arbitrary in a classical 
theory. 
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uniformly bounded, < 1 — | !>[/□< 2, due to 
the ST/ (2) compactification. Thinking in terms of 
statistical mechanics for our classical lattice sys- 
tem, we expect that after some time the typical 
field configuration has equally much energy in all 
modes of vibration - independent of the frequency| 
- and the total amplitude of the classical field, and 
the energy per lattice plaquette, is thus small for a 
fixed low energy. For low energy, when the average 
energy per plaquette is small, the lattice artifacts 
due to the compactness of the gauge group are 
thus negligible. 

(2) For small energy per plaquette, we thus 
expect that the dominant form for lattice artifacts 
is due to the fact that an appreciable amount of 
the activity (for example the energy) is in the field 
modes with wavelengths comparable to the lattice 
constant a. This short wavelength activity at lat- 
tice cut-off scales is unavoidable in the limit of 
long time simulation of an initially smooth field 
configuration (relative to the lattice spacing), or 
already after a short time if we initially have an 
irregular field configuration. 

In which way can the classical lattice gauge 
theory approach a continuum limit ? This is a 
difficult question to which we shall only be able to 
provide a very partial answer here. 

For the study of the time evolution of Yang- 
Mills fields which initially are far from an equi- 
librium situation, the (classical) field modes will 
exhibit a never ending dynamical cascade towards 
the ultraviolet and after a certain transient time, 
the cut-off provided by the spatial lattice will pre- 
vent the lattice gauge theory from simulating this 
cascade. It is therefore immediately clear that 
the lattice regularized, classical fields will not ap- 
proach a "continuum limit" in the sense of simulat- 
ing the dynamical behavior of the classical contin- 
uum fields in the t — ► oo limit. For the simulation 
of classical continuum gauge fields far from equi- 
librium, we conclude that we will have the best 
"continuum limit" if we simulate, for a short pe- 
riod of time, an initial smooth ansatz[] for the fields 

e Note, this situation is very different in the quantum 
case. Planck's constant h introduces a relation E = hui be- 
tween the energy of a mode of vibration and its frequency, 
implying that a mode with a high frequency also has a high 
energy. With a given available finite total energy, modes 
with high frequencies will therefore be suppressed. Quan- 
tum mechanically, we thus have that at low energy only 
excitations of the longest wavelengths appear. 

^By "smooth" configurations we mean lattice config- 



in the region of low energy per plaquette|| 

Gauge fields exist, however, as quantum the- 
ories in Nature and the interesting definition - as 
concerns applications in physics - of a "continuum 
limit" of the classical lattice gauge theory, is thus 
to identify regions in the parameter space for the 
classical lattice gauge theory which probe the be- 
havior of the time evolution of semi-classical initial 
configurations (with many quanta) of the quantum 
theory for a shorter or longer interval of time or in 
an equilibrium situation. (See also e.g. discussion 
in RefEa). Since the classical lattice gauge the- 
ory does not contain a relationship like E = Hlu 
(implying a damping of the high frequency modes 
relative to the soft modes), we must expect that 
the simulation of quantum gauge theory will be 
distorted by this fact (even if implemented with 
effective lattice Hamiltonians as e.g. devised by 
RefET). 

We shall in the following restrict attention to 
real-time simulations of the classical lattice gauge 
theory which have reached an equilibrium situa- 
tion on the lattice. Such simulations are hoped 
u to yield insights into the dynamical behavior of 
(long wavelength modes in) the high temperature 
Q.C.D. fields in situations where an equilibrium 
situation has been reached. 

Before we attempt an analysis of aspects of 
chaos (in time and space) of the lattice gauge the- 
ory let us first note that it is far from obvious how 
to establish contact between an effective lattice 
temperature of the classical lattice gauge theory 
(which has reached an equilibrium situation due 
to the presence of the lattice cutoff A = 1/a) and 
the physical temperature T of the Q.C.D. theory 
which is a quantum theory and which can reach an 
equilibrium situation due to a cut off of quantum 



urations which are 'good' approximations to continuum 
configurations. 

9 I.e. we are here imagining a situation where the spa- 
tial correlations in the monitored field variables are so large 
that they lose memory of the underlying lattice structure 
(including the lattice spacing a). In the extreme opposite 
limit, one could imagine situations with randomly fluctu- 
ating fields on the scale of the lattice constant, i.e. with 
(almost) no spatial correlations from link to link. If the 
field variables fluctuate independently of each other (inde- 
pendent of their neighbors), one could imagine the model 
to be invariant (with respect to the monitoring of many 
variables) under changes of the lattice spacing, a. Thus, 
it appears that lattice cut-off independence of numerical 
results can not be a sufficient criterion for the results to 
report "continuum physics". 
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mechanical origin (i.e. ultimately due to the exis- 
tence of a Planck xxmstant K). Some discrepancy 
in the literature E2lE3n3 illustrates that this is not 
an easy question. 

In the classical lattice gauge theory we have 
an "effective lattice temperature" Tl = 1/(3 where 
(3 has been determined by looking at the probabil- 
ity distribution (Gibbs distribution) of the kinetic 
energy (or the magnetic energy in a plaquette), 
p{EkX oc exp(— (3Ek) (after calibration) . Cf. e.g. 
Ref. I (p. W>\^) and Ref. M. 

In RefBila'Ea it is asserted that the physical 
temperature T (of the quantum Q.C.D. fields in an 
equilibrium situation) and the average energy En 
per plaquette in the classical lattice gauge study 
in equilibrium are related by 



Ec 



(n 2 -l)T (for SU(n)). 



(5) 



It is not clear to us how serious one should take 
this relationship (one objection being that a scale 
has not been fixed in the continuum limit of the 
classical theory) but we shall leave a discussion of 
this issue aside in the following. 

Let us now attempt an understanding of 
chaotic aspects of the lattice gauge theory - i.e. 
the temporal chaos (as monitored by the spectrum 
of Lyapunov exponents) and its relationship with 
spatial chaos (as monitored by the spatial corre- 
lations, e.g. a spatial correlation length) in the 
dynamical system. ___ 

A sequence of articles^Olia and a recent book 
by Biro, Matinyan and Mullefl present numerical 
results for the classical SU{2) lattice gauge model 
and provide evidence that the maximal Lyapunov 
exponent is a monotonically increasing continuous 
function of the scale free energy/plaquette with 
the valufipZpxD at zero energy. 

RefErlLj'Ej reports a particularly interesting in- 
terpretation of numerical results for the dynamics 
on the lattice, namely that there is a linear scaling 
relation between the scale free maximal Lyapunov 
exponent, X m ax(a = 1) and the average energy per 
plaquette E n (a — 1). The possible physical relet 
vance of this result follows from the observation!^ 
that when we rescale back to a variable lattice 
spacing a we note that the observed relationship 
is in fact a graph of a\ max {a) as a function of 
a E a (a). Thus being linear, 



cancellation of a factor a implies that X max (a) = 
const x En (a) and thus there is a continuum limit 
a — > 0, either of both sides simultaneously or of 
none of them. In the particular case where the 
energy per mode (oc energy per plaquette) is in- 
terpreted to be a fixed temperature T (cf. equa- 
tion © above), one deduces that the maximal 
Lyapunov exponent has a continuum limit in real 
time, proportional to the temperature of the gauge 

Note, also, the suggested relationshi]S'EE0 be- 
tween the maximal Lyapunov exponent A of the 
gauge fields on the lattice and the "gluon damp- 
ing rate" 7(0) for a thermal gluon at rest, arrived 
at in re-summed perturbation theory in finite.te.rrij- 
perature quantum field theory, see also RefsESEa. 
For the SU(2) gauge theory, this suggested rela- 
tion reads 

A = 2 7 (0) = 2 x 6.635-^- 5 2 T - 0.34 g 2 T (7) 

247T 

It is also our understanding-.that 7(0) is a quan- 
tity of semiclassical originliJ A relation like (Q) 
is nevertheless remarkable in chaos theory since 
it suggests that a complicated dynamical quan- 
tity like a temporal Lyapunov exponent (which is 
usually only possible to extract after a consider- 
able numerical effort) is analytically calculable by 
summing up some diagrams in finite temperature 
quantum field theory. 

In Refta we argued (in-pfissL of the numerical 
evidence presented in e.g. Stalls) that the appar- 
ent linear scaling relation (|^) is a transient phe- 
nomena residing in a region extending at most 
a decade between two scaling regions, namely 
for small energies where the Lyapunov exponent 
scales with an exponent which could be close to 
k 1/4, and a high energy region where the scal- 
ing exponent is at most zero: 



dlog A 
dlogE 



= fc~l/4 for£^0 
< for E -> 00 



(8) 



a Xr 



const x a E a (a) 



(6) 



The proposed scaling relation k ~ 1/4 was, for 
reasons we shall give below, suggested to hold in 
ic limit E — > from general scaling arguments 
of the continuum classical Yang-Mills equations 
in accordance with simulations on homo gen eous 

aodels and consistent with the figures inEj^md 
. However, a more recent numerical analysis^ ar- 
gues rather conyjrLcingly that the data points pre- 
sented in e.g. mm were subject to "finite time" 
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artifacts, and that long time simulations with a 
correct procedure for extracting the principal Lya- 
punov exponent could support that k = 1 even in 
the limit as E — » 0. 

Before we discuss the limit of low energies 
(E — > 0) let us note that regarding the limit for 
high energieager plaquette a rigorous resultQ (also 
reported in no) shows that the SU (2) scale free 
(a = 1) lattice Hamiltonian in d spatial dimen- 
sions has an upper bound for the maximal Lya- 
punov exponent 

X B = yj(d- l)(4 + \/l7) (9) 

which for d = 3 becomes = 4.03... This result 
is arrived at by constructing an appropriate norm 
on the phase space and showing that the time 
derivative of this norm can be bounded by a con- 
stant times the norm itself, hence giving us an up- 
per bound as to how exponentially fast the norm 
can grow in time. The upper bound (^) shows that 
a linear scaling region, i.e. rj = dlog X/d\ogE = 1, 
cannot extend further than around E ~ 10 on 
the figure 1. Beyond that point the maximal Lya- 
punov exponent either saturates and scales with 
energy with an exponent which approaches zero 
or it may even decrease over a region of high ener- 
gies, yielding a negative exponent thus justifying 
the E — ► oo limit of equation (g). 

It should be noted that the upper bound (|) is 
independent of the lattice size and the energy (but 
scales with 1 /a) . There is quite a simple intuitive 
explanation for the saturation of the maximal Lya- 
punov exponent in the regime for high energy per 
plaquette of the SU (2) lattice gauge model, since 
the potential (magnetic) term in the Hamiltonian 
(|) is uniformly bounded, < 1 - \TrU a < 2, due 
to the SU(2) compactification. For high energies, 
almost all the energy is thus put in the (integrable) 
kinetic energy term in the lattice Hamiltonian (|]), 
and the spectrum of Lyapunov exponents will sat- 
urate as the energy per plaquette increases for the 
model. (The chaos generated by the non-linear 
potential energy term does not increase, only the 
energy in the kinetic energy (the electric fields) 
increases) . 

As we sh all argue (a discussion which is rather 
suppressed inEjO'ESp), it is however in the oppo- 
site limit, i.e. in the limit where the average energy 
per plaquette goes to zero, i.e. E — ► 0, that the 
equilibrium lattice gauge theory has the possibil- 



ity to reach contact with "continuum physics" - 
not in the sense (as we have seen) that the lattice 
system simulates the original continuum classical 
gauge fields (over a long time span), but in the 
sense that the lattice gauge theory - considered as 
a statistical system in its own right - will develop 
spatial correlations in the fields (as monitored by 
a correlation length £, say) which may approach a 
formal "continuum limit" , 

f -> oo as E -> (10) 

There is a well known analogjEl between a quan- 
tum field theory in the Euclidean formulation with 
a compactified (periodic) imaginary time axis < 
t < and finite temperature statistical mechan- 
ics of the quantum field theory at a temperature 
T = 1/(3 which is inversely proportional to the 
above-mentioned time extension. 

Calculations in finite temperature Euclidean 
quantum field theory suggest that a characteris- 
tic correlation length of static magnetic fields in 
the thermal quantum gauge theory is of the order 
£ ~ (g 2 T) _1 . If we (cf. discussion above, equa- 
tion (j|)) make use of the relationship, E a = 2T 
(for SU(2) Yang-Mills theory) between tempera- 
ture T and average energy per plaquette E a , this 
could suggest that the (magnetic sector) of field 
configurations on the lattice will have characteris- 
tic correlation lengths of the order 

£ ~ (.g'T)- 1 ~ {-g 2 E n )- X -> oo as E a -» . 

(11) 

In systems with space-time chaos (and local prop- 
agation of disturbances with a fixed speed ~ c) 
one often has a relationship between a correlation 
length (coherence length) £ and the maximal Lya- 
punov exponent A (see also e.g. discussions in 
RefO) of the form £ ~ c/A. This suggests, via 
equation (pr|), 

A ~ c/£ ~ g 2 E n (12) 

i.e. a Lyapunov exponent which scales linearly 
with the average energy per plaquette E a . Clearly, 
it deserves further investigation to establish a pre- 
cise relationship between X ma x and the spatial cor- 
relation length £ (in lattice units). However, we 
emphasize that it is in the limit as the energy per 
plaquette goes to zero (and where the lattice Lya- 
punov exponent X m ax goes to zero) that we expect 
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the spatial correlation length £ (in lattice units) 
will diverge. The finite size of the lattice makes it 
difficult to analyze the behavior of the gauge fields 
and the principal Lyapunov exponent on the lat- 
tice in the limit E — > 0. 

Let us reiterate the arguments which led us to 
the conclusion (in Refta) that the scaling exponent 
of the Lyapunov exponent likely will be closer to ~ 
1 /4 in the limit E — > than the scaling exponent 
~ 1 observed in the intermediate energy region 
1/2 <E <4. 

In figure 1 we displav^a log-log plot of the re- 
sults obtained in RefBE3il3. It appears that for 
E ~ 1/2 there is a cross-over to another scal- 
ing region. Although the data points in this re- 
gion are determined with some numerical uncer- 
tainty (finite time artifacts, as argued in Rem) we 
did note that they are consistent with a scaling 
with exponent ~ 1/4. Moreover, it is well known 
that the homogeneous Yang-Mills equations have a 
non-zero Lyapunov exponent which by elementary 
scaling arguments scales with the fourth root of 
jtlje energy density and has the approximate form 



A„ 



1/3 £ 1/4 



(13) 



On the lattice the above scaling relation is valid 
for spatially homogeneous fields, i.e. the max- 
imal Lyapunov exponent scales with the fourth 
root of the energy per plaquette. By continu- 
ity, fields which are almost homogeneous on 
the lattice will, in their transient, initial dynami- 
cal behavior, exhibit a scaling exponent close to 
1/4. In fact, the same scaling exponent would 
also hold for the inhomogeneous Yang-Mills equa- 
tions U had the fields been smooth relative to the 
lattice scale, so derivatives, d^, are well approx- 
imated by their lattice equivalent and we are al- 
lowed to scale lengths as well. This is seen from 
scaling arguments for the continuum Yang-Mills 
equations D^F^ — where = <9 M — iglA^, \ 
and F" v = d^A u -d u A^-ig[A",A u ]. Not taking 
boundary conditions into account, these equations 
are invariant when and A^ are scaled with the 
same factor a. That is, if A(x,t) is a solution to 
the equations, then ^A(ax, at) is also a solution. 
The energy density E, which is quadratic in the 
Yang-Mills field curvature tensor F^, then scales 



h Such field configurations are (ungeneric) examples of 
field configurations which exhibit correlation lengths much 
larger than the lattice spacing. 



with of 4 . This indicates @ that if we perform a 
measurement of the maximal Lyapunov exponent 
A mol over a time short enough for the solutions to 
stay smooth, then \ max scales with E 1 / 4 . Such a 
scaling was also observed for "smooth" configura- 
tions by Miiller and Trayanov p. 3389 in RefO. 
These scaling arguments do not carry over to in- 
finite time averages since solutions on the lattice 
tend to be irregular. 

In any case we expect that the lattice gauge 
theory does not probe field configurations which 
are substantially correlated beyond one lattice 
unit when the lattice gauge theory is monitored for 
energies per plaquette in the regime E ~ 1/2 — 4. 
That is, we expect £/a ~ 1 for E ~ 1/2 — 4 
and thus the system is indeed very far from reach- 
ing contact with "continuum physics" when mon- 
itored in that regime of energies. For example, 
the numerical results for the Lyapunov exponents 
as presented e.g. in Ref. BEjIia (c.f. e.g. fig. 
8.4 in Rem) report on lattice size N = 20 for en- 
ergies E ~ 1/2 — 4 and coupling strength g w 2, 
which suggest (cf. also the equation (|ll|)) that the 
lattice gauge theory in this regime of parameters 
monitor lattice field configurations which have a 
correlation length up to a few lattice units only. 
There is also numerical evidence | that the corre- 
lation length for energies E ~ 1 is of the order of a 
few lattice units. This observation applies also to 
the numerical studias-.of the SU(3) lattice gauge 
theory reported inlilM in the range of energies per 
plaquette E a ^4 — 6. 

These observations are further substantiated 
by the numerical simulations of Reflia (figure 12) 
for the lattice gauge theory with a U(l) group 
showing a steep increase of the maximal Lya- 
punov exponent with energy/plaquette in the in- 
terval 1 < E < 4. The continuum theory here 
corresponds to the classical electromagnetic fields 
which have no self-interaction, and thus the Lya- 
punov exponent in this limit should vanish. .-The 
discrepancies in this case were attributed inEa to 
a combined effect of the discreteness of the lattice 
and the compactness of the gauge group U(l) and 
were not connected with finite size effects. This 
suggests strongly that we, in the case of numer- 
ical studies of a SU(2) Kogut-Susskind Hamilto- 

1 J. Ambj0rn, priv. communication. See also e.g. Ref. 
which reports investigations of finite temperature QCD on 
the lattice with measurements of a magnetic mass (inverse 
correlation length) of the order m mag ~ 0.5 g 2 T. 
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Figure 1: The maximal Lyapunov exponent as a function 
of the average energy per plaquette for the SU (2) and the 
U(l) lattice gauge theory. The data points (diamoifdp for 
SU(2), tictosifor U(l)) are adapted from Miiller et al. £3 and 
Biro et al. p. 192. The solid line is a linear fit through 
the origin, the dashed line is the function ^ X 1/3 X i? 1 / 4 
(half of the homogeneous case result (H)). The dot-dashed 
line shows the rigorous upper bound for the SU (2) lattice 
model (saturation of temporal chaos). As is seen, the linear 
scaling region for the SU(2) data is positioned where the 
(7(1) data display strong lattice artifacts. 

man system, cannot base continuum physics on 
results from simulations in the same interval of 
energies where the U(l) simulations fail to display 
continuum physics. On the contrary, we suspect 
that what could reasonably be called "continuum 
physics" (for the lattice gauge theory considered as 
a statistical system in its own right) has to be ex- 
tracted from investigations of the Kogut-Susskind 
lattice simulations for energies per plaquette which 
are at least much smaller than E ~ 1/2. 

We would like to conclude with some final 
points of discussion. 

Exactly in which way we may speak about 
a "continuum limit" in a simulation of a classi- 
cal gauge theory in real time on a spatial lat- 
tice (a theory which does not have a continuum 
limit without lattice regularization) appears to be 
a question of some fundamental as well as practi- 
cal interest, since not many non-perturbative tools 
are available to study the time evolution of quan- 
tized Yang-Mills fields. 

Sio£c by intrinsic scaling arguments (cf. 
Refs, uu) one has a functional relation between 



aA maa; (a) and aE a {a), a scaling (at small energies 
E — » 0) with exponent k, according to the scaling 
relation (||), would imply that 

A(a) cx a k - 1 E u {a) k , (14) 

in which case one cannot achieve a continuum limit 
simultaneously for the maximal Lyapunov expo- 
nent and the temperature (assuming that it is 
proportional to the average energy per plaquette 
E(a), as given by equation (||)), except in the spe- 
cial case where k = 1. In particular, if k < 1, the 
former would be divergent if the temperature is 
kept fixed. There is no particular contradiction in 
this statement, however, as there is - a priori - no 
reason for having a finite Lyapunov exponent in 
the continuum limit. The erratic and fluctuating 
behavior of the fields one expects in time as.. well 
as in space (for numerical evidence, cf. alsoEJ) on 
very small scales could suggest that a Lyapunov 
exponent would not be well defined in the "con- 
tinuum limit" (as a — > and E — > 0). Clearly, this 
question deserves further investigation. 

If the scaling-relation k = 1 according to equa- 
tion (||) holds in the limit E — > for the lattice 
gauge theory, it will in a most striking way illus- 
trate the point that the continuum gauge theory 
(with scaling k ~ 1/4) and the lattice gauge the- 
ory (probed in a situation where it has reached 
equilibrium) are two entirely different dynamical 
systems - despite the lattice theory (^) is at first 
set up to be an approximation to the continuum 
theory. (Thus the lattice theory does not simu- 
late the continuum theory; it is an entirely new 
statistical theory in its own right and - important 
for applications in physics - the relationship with 
quantum non-Abelian gauge theory remains to be 
established on a more rigorous basis). Under- 
standing of this crucial difference could, perhaps, 
be obtained from renormalization group analysis: 
We may say that the limit E — » is a critical point 
for the lattice Hamiltonian (|]), i.e. that the lat- 
tice correlation length diverges in that limit. As 
is common for field theories studied in a neigh- 
borhood of a critical point, one expects classi- 
cal scaling arguments to break down, or rather, 
that scaling relations are subjected to renormal- 
ization which gives rise to anomalies in the scaling 
exponents. Often such anomalous scaling expo- 
nents seem to be 'ugly' irrational numbers, per- 
haps with the 2-d Ising model as a notable ex- 
ception. It would therefore be quite miraculous 
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if the 1 /4 classical-scaling of the Lyapunov expo- 
nent emphasized above renormalizes to an expo- 
nent which equals unity. Even if this did happen 
it is not clear that such a behavior should be inde- 
pendent of the regularization procedure employed, 
since Lyapunov exponents are measures of local in- 
stabilities, i.e. short range rather than long range 
structures. 

As regards the numerical evidence for the re- 
lationship (^) between the maximal Lyacunov ex- 
ponent and the "gluon damping rate" 23 we note 
that in the hierarchy of scales g 2 T <C gT <C T 
(this separation of scales is assumed in hot per- 
turbative gauge theory) the "gluon damping rate" 
is connected to the decay of the "soft" modes 
~ g 2 T. In order for the lattice gauge theory to 
probe decays of "soft" gauge modes, this requires 
the existence of some "soft" modes on the lat- 
tice (in a background of "hard" modes). Thus 
we must monitor the lattice gauge system in a re- 
gion where there are fields (in equilibrium with 
"hard" modes A ~ 1/a) which have spatial corre- 
lations substantially larger than the lattice unit 
a. As we have seen this is not the case in the 
regime where the lattice gauge theory was, studied 
(N = 20, g w 2,E ~ 1 - 4) in Ref. (Mfca pro- 
viding the numerical support for equation (^). It 
appears, that the numerical support of the relation 
([?]) is somewhat of an accident if the relation ([?]) is 
to be interpreted as a "continuum result" . Studies 
of the lattice gauge theory have to be conducted 
for much smaller energies per plaquette which are, 
however, also difficult, since finite size effects then 
will become of appreciable size. •[] 
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